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The detailed origin of the diffuse gamma-ray background is still unknown. However, the contribution of unresolved 
sources is expected to induce small-scale anisotropies in this emission, which may provide a way to identify and 
constrain the properties of its contributors. Recent studies have predicted the contributions to the angular power 
spectrum (APS) from extragalactic and galactic dark matter (DM) annihilation or decay. The Fermi-LAT collabo- 
ration reported detection of angula r power with a significance larger than 3cr in the energy range from 1 GeV to 10 
GeV on 22 months o f data [Ackermann et ah 2012 . For these preliminary results the already published Fermi-LAT 
APS measurements Ackerma nn et al.||2012| are compa red to the accurate predictions for DM anisotropies from 
state-of-the-art cosmological simulations as presented in [Fornasa et al.|2013] to derive constraints on different DM 
candidates. 



1. THE DIFFUSE GAMMA-RAY 
BACKGROUND 

The diffuse gamma-ray background is characterized 
by an isotropic or nearly isotropic distribution and 
is therefore known as the Isotropic Gamma-ray Back- 
ground (IGRB) Kraushaar et al.|1972 . It is constituted 
by gamma rays produced by various sources, including 



blazars, pulsars, and possible DM structures, not yet de- 
tected due to the limited angular resolution and photon 
statistics of the Fermi-LAT. Figure [l] shows the IGRB 
spectrum and the estimated contributions from unre- 
solved blazars, star-forming and radio galaxies. The an- 
gular distribution of photons in the diffuse background 
may contain information about the presence and the 
nature of these unresolved source populations. 
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Figure 1: IGRB spectrum and the contribution from all the different extragalactic source classes that have been detected 
by Fermi. The golden band shows the sum of all the source populations. By summing all contributio ns it is appa rent that 
there is still room for other components at all energies within the uncertainties Sanc hez-Conde||2012| . See also [Stecker 
and Venters||2011| for other theoretical estimates of the relative contributions of unresolved blazars and star-forming 



galaxies to the IGRB. 



2. FERMI-LAT MEASUREMENT OF THE 
ANGULAR POWER SPECTRUM IN THE 
IGRB 



In |Ackermann et aL 2Q12| the first 22 months of 
Fermi-LAT data were analyzed, dividing the energy 
range between 1 GeV and 50 GeV in 4 energy bins. 
The point sources in the first year catalogue Abdo et al. 
2Q1Q| have been masked, as well as the emission within 



a band of 30 degrees above and below the Galactic plane 
. The masking was done to cover the regions in the sky 
where the emission is dominated by resolved sources and 
by the Galactic foreground, and to restrict the analysis 
only to where the IGRB is a significant component. 



Two definitions for anisotropics were used in Acker- 
mann et al.]|2012| : 



• Intensity APS: An intensity map can be de- 
composed in spherical harmonics, 



(1) 



where coefficients ai m determine the APS which is 
given by C\ = ^|a/ m | 2 ^. This definition is partic- 
ularly useful because it gives us the dimensionful 
size of intensity fluctuations and can be compared 
with predictions for source classes whose collective 
intensity is known or assumed. 

• Fluctuation APS: can be derived from the in- 
tensity APS, dividing by the average inten- 
sity squared. The fluctuation APS is energy- 
independent for a single source class, if all mem- 



bers of the source class share the same observed 
energy spectrum. 

The Fermi-LAT collaboration reported detection of 
angular power in all 4 energy bins considered, with a 
signicance larger than 3a in the energy bins from 1 GeV 
to 10 GeV. The data have been compared with the APS 
of a source model made of i) the point sources in [Abdo 
et al.|2010 , ii) a model for the interstellar diffuse emis- 
sion and iii) an isotropic co mponent at the level of the 
IGRB in Abdo et al. |2010|. The model angular power 



at 155 < I < 504 is consistently below that measured in 
the data. 

Despite the mask applied along the Galactic plane, 
some known Galactic emission can extend to high lat- 
itudes. Therefore a model of the Galactic foregrounds 
was subtracted from the data, and then the APS of the 
residual maps was calculated. This measure ment is re- 
ferred to as the cleaned data in Ackermann et al.|2012 



We use this second measurement in this work. 



3. DARK MATTER PREDICTIONS 

The APS of gamma rays from DM annihilations or 
decays has been c omputed from the al l-sky template 

The authors 



maps produced in Fornasa et al. 



2013 



of [Fornasa et al.||2013| used the Millennium-II N-body 
simulation to model the abundance and the clustering 
of extragalactic DM halos and su bhalos. The technique 
presented in [Zavala et al.||2010 , based on the random 
repetition of copies of the Millennium-II simulation box, 
is implemented to probe the universe up to z = 2. The 
emission from DM halos with a mass below the resolu- 
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tion of the simulation was estimated assuming that the 
halo number density and the mass-luminosity relation 
obtained from the halos in Millennium-II remains un- 
changed below the resolution, down to the minimal self- 
bound halo mass M rain . On the other hand, the contri- 
bution of low mass sub halos was modeled following the 
technique described in 



Kamionkowski et al. 
Sanchez- Conde et al.| |2011 



2010 



and 

extended in 

The smooth DM halo of the Milky Way was 
parametrized in Fornasa et al.|2013 as an Einasto pro- 
file, since this provides the best fit to the Milky Way- 
like halo obtained in the Aquarius N-body simulation. 
Galactic subhalos down to 10 5 M were accounted for 
directly from the Aquarius simulation, while we use the 
same procedure as before to account for the contribu- 
tion of unresolved subhalos. It has been shown that 
such objects do not contribute significantly to the total 
intensity APS. 

In |Fornasa et al.||2013 the effect of the assumptions 
made in the modeling of the DM distribution also was 
estimated, looking for their effect both on the intensity 
of the DM-induced emission and on its APS. The two 
most relevant sources of uncertainty are the amount 
of substructures hosted by DM halos and the value of 
Mmin. The first gives rise to an uncertainty of a factor 
20-30 both in the average intensity and in the intensity 
APS; when the uncertainty on the value of M m ^ n is 
taken into account, these factors build up to 40 and 
100 for the average intensity and the intensity APS, 
respectively. 



4. SETTING CONSTRAINTS 



Here we present the method used to set conservative 
limits on the thermally averaged cross section, (cry), 
for DM annihilation into three different channels, bb 
quarks, and r + r~ leptons. 

• To set constraints we use the foreground- cleaned 



C p shown in table II of Acker mann et al. 2012 



There are four C p values measured corresponding 
to four energy bins, 1 - 2 GeV, 2-5 GeV, 5-10 
GeV, and 10 — 50 GeV. We use them indepen- 
dently to set limits. 



• We know that the IGRB anisotropy has multiple 
contributors, therefore these constraints are con- 
servative. Other contributors to IGRB anisotropy 
are not well known, but we already have con- 
straints on the contribution of blazars Cuoco e1 
al.||2Q12 . We subtract this contribution from the 
measured C p and require that the DM-induced 
APS does not overshoot this new limit. 



5. PRELIMINARY RESULTS 

Figures [2] -[4] show the 95% CL limits on the annihila- 
tion cross section for three different channels, bb quarks, 
and r + r~ leptons, respectively. The main un- 



certainty in the predictions obtained in |Fornasa et al.| 
|2013| lies in the properties of low- mass subhalos, be- 
low the mass resolution of the simulations. Different 
values of the "subhalo boost" strongly affect the pre- 
diction for the DM-induced gamma-ray intensity and 
its anisotropics. Fornasa et al. 2013 considered two 



benchmark scenarios for subhalos, assuming that the 
uncertainties can be modeled by changing the subhalo 
abundance: i) the LOW case, where halos are rela- 
tiv ely poor in subhalos, acco rding t o the predictions 

2010| and [Sanchez- Conde et 



of Kamionkowski et al. 



al. |2011 ; the constraints using this scenario are shown 
in the left-upper panel of figures [2] -|4j and ii) the HIGH 
case, with large subhalo boosts, compatible with wha t 
was found by [ Pinzke et aL]|2011[ |Gao et al.||2012|; 
the right-upper and lower panels in the figures [2] - |4| 
use this case. In the lower panels the predictions for 
the DM-induced APS (with a HIGH subhalo boost) are 
compared with the APS measured by Fermi-LAT once 
the contribution from blazars is subtracted. This rep- 
resents the case where the most stringent constraints 
on (av) are obtained. We are currently updating the 
Fermi-LAT APS measurement using ~ 45 months of 
data, we expect to improve the constraints presented in 
these proceedings after include the new measurement in 
the method shown in section 3J 



• Requiring that the DM-induced intensity APS av- 
eraged in 155 < I < 204 does not overshoot the 
measured C p in the 155 < I < 504 multipole range 
plus 1.64 times its error leads to 95% CL limits 
on (0-v) 1 . 



1 We assume Gaussian erros, then the value of 1.64 is based on 
the fact that 95% of the area of a Gaussian distribution is within 
1.64 standard deviations of the mean. 
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Figure 2: 95% CL limits on the annihilation cross section for the bb channel. 
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Figure 3: 95% CL limits on the annihilation cross section for the channel. 
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Figure 4: 95% CL limits on the annihilation cross section for the r + r channel. 
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